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We investigate the effect of weak randomness on the antiferromagnetic anisotropic spin-1 chain. 
We use Abelian bosonization to construct the low-energy effective theory. A renormalization group 
calculation up to second order in the strength of the disorder is performed on this effective theory. 
We observe in this framework the destruction of the antiferromagnetic ordered phase a la Imry-Ma. 
We predict the effects of a random magnetic field along z axis, a random field in the XY plane as 
well as random exchange with and without XY symmetry. Instabilities of massless phases appear 
in general by mechanisms different from the case of the 2-leg spin ladder. 
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I. INTRODUCTION 

Quantum one-dimensional spin systems are characterized by a variety of behavior much richer than their higher- 
dimensional counterpart. For example, antiferromagnetic spin chains have different physical properties according to 
whether the spins are integer or half-integerEI. Let us consider the prototypical Heisenberg spin chain : 

W = J^Si-S i+1 , (1.1) 

i 

where i are the sites along the chain and Si quantum spin operators with Sj 2 = S(S + 1). The exchange coupling J 
is positive, i.e. antiferromagnetic. Then the S=l/2 chain is gapless and has algebraic decay of the spin correlations 
(S • S„) ~ {-) n /n for n ->■ oo. 

The S=l chain has a quite different physics : the isotropic antiferromagnetic chain is in the Haldane phase. There 
is a gap for spin excitations and the correlations decay exponentially. When adding anisotropy, this Haldane phase 
survives in a region of parameters around the isotropic pointnra. When anisotropics become too large, then the system 
is eventually driven towards more conventional ferro or antiferro phaseso. 

This spin-1 chain is closely related to the two-leg spin-1/2 ladder. In the ladder system, one is dealing with 
two S=l/2 spin chains coupled by rungs perpendicular to the chains. The ladder has a gap for all values of the 
perpendicular coupling. In fact, it is known that the S=l chain and the two-leg spin ladder share the same low-energy 
effective field theorytffl. 

Understanding the influence of disorder on such spin systems is an important issue. In all real materials there is 
of course some amount of disorder. One of the probe of the physics of the system is to introduce artificially defects 
or impurities. The study of disorder in one-dimensional quantum systems is also of considerable theoretical interesifl. 
Recently it has been realized that it is possible to obtain exact results in such disordered systems, while this is 
generally impossible in the classical realm. 

In the case of the spin-l/2jchain, one simple kind of disorder is a distribution of quenched bond disorder. The real- 
space renormalization groupQ study shows that the ground state is a collection of singlet bonds over arbitrarily large 
distances in a random pattern according to the initial distribution of randomness. This is the so-called random singlet 
phaseo. The influence of weak disorder on the spin-1/2 chain has been also studied by means of a renormalization 
group (RG) study on the effective field theorya. This method is a perturbative RG calculation but it is easy to study 
widely different kinds of disorder. The real-space RG method has .also, been applied to systems with a distribution of 
ferro- and antiferromagnetic bondsEll as well as dimerized systemsExt 2 ! There are also direct numerical studiestia 

In the case of the spin-1 chain we are dealing with a gapped system that has alsOpSome hidden long-range order 
as exhibited by the approximate valence-bond-solid (VBS) ground state wavefunctionsotl Certainly one expects the 
influence of disorder to be quite differentlia. A real-space RG has been performed on the isotropic S=l chairJl3. It 
shows that there is destruction of the Haldane gap beyond a nonzero disorder strength, then there is a gapless phase 
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which sustains the hidden long-range order and then for even stronger disorder one finds the featureless random-singlet 
phase. This is analogous to the doping of conventional superconductors with magnetic impuritieslHI : here the hidden 
long-range order plays the role of the superconducting condensate. 

In this paper we study the S=l spin chain with exchange and single-ion anisotropics. This model has a rich phase 
diagram and the influence of disorder is extremely sensitive to the nature of the ground state of the unperturbed 
system. We use a perturbative RG treatment of the disorder on the effective field theory. The average over the 
disorder is performed by use of the replica trick as pioneered by Giamarchi and SchulzGJj in the context of the one- 
dimensional electron gas. Here we perform the-calculation up to second order in the disorder strength. This is 
important as shown by Fujimoto and Kawakamiij because, at this order of perturbation, one captures the interplay 
between interactions and disorder when they are both relevant. In the spin- 1/2 case they have shown that it is 
possible to explicitly obtain from the RG treatment the destruction of the antiferromagnetic phase under a random 
z-axis magnetic field for an .arbitrarily small amount of randomness as predicted by the argument of Imry and MaEJ. 
The case of the spin ladderEil is closely related to our work but the relevant operators are not always the same. In 
the spin-1 chain there are operators generated by disorder that involve products of the two spin- 1/2 operators used 
to describe the spin-1 degrees of freedom. They do not exist in the ladder case. 

In section II, we give the construction of the effective bosonic theory for anisotropic S=l chain without randomness. 
In section III, we translate in the bosonic language, the operators induced by the disorder. Section IV shows why 
gapped phases are in general stable under weak disorder. Section V contains the derivation of the renormalization 
group equations. Section VI gives the results of our study on each phase of the S=l chain. Finally section VII contains 
our conclusions. 



II. THE PURE S=l SPIN CHAIN 



In this section we give the construction of the effective Bosonic low-energy theory for the S=l spin chain. We will 
use this continuum theory as a starting point for the treatment of randomness in sect. III. We concentrate on the 
anisotropic Heisenberg S=l antiferromagnetic spin chain with the following Hamiltonian : 

H = J2 S ? S ?+i + S?S? +1 + AS*Sf +1 + D{S*f. (2.1) 

i 

The Sf are S=l quantum spin operators, A is the exchange anisotropy and D is the single-ion anisotropy_ All energies 
are measured in units of J. The phase diagram of this model is well-known in the absence of disorder!]. There are 
two phases with conventional long-range magnetic order : a ferromagnetic (F) phase for A large and negative and 
an antiferromagnetic (AF) phase for A large and positive. In the intermediate regime, there is the Haldane phase 
(n) which is gapped and has no obvious long-range order. In particular it does not break translational or rotational 
symmetry. There are two gapless XY-like phases that differ only in the behaviour of the spin correlations. Finally, 
when D is very large and positive there is a non-magnetic phase which is trivially obtained by perturbation theory 
from the strong-coupling limit where the ground state has all spins set to Sf = : 1 000 . . .00), called the large-D 
phase (D) in what follows. 

The starting point is to write each spin-1 operator as the sum of two spin-1/2 operators Sf = erf ^t", an 
approximation suggested by Luther and PescheEj. It has been shown by Schulz that Abelian bosonizationo is able 
to reproduce most of the phase diagram. One first performs a Jordan- Wigner transformation on each spin-1/2 a and 
t by introducing two species of spinless fermions. Explicitly we write : 

( _)n c t e -£ m <»<4c m) (2 . 2a) 

(-)"c„e"£ m <„^" ; (2.2b) 

4cn - \ , (2.2c) 

where c is a fcrmion operator and similarly one introduces a d fermion for the other set of spins-1/2 r". Since 
one is interested in the long- wavelength, low-energy behaviour, it is then convenient to take the continuum limit 
by introducing two continuous Fermi fields ipc( x ) and ipd(x)- The Hamiltonian ( [2.l| ) leads then to a theory of two 
interacting Fermi fields in one space dimension. The physical content of this Fermi system is obtained by translation 
into the language of interacting BosonscJ. In fact, in one space dimension a theory involving a Fermi field ip( x ) can 
be translated in a Bose field theory by the following relations : 
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(2.3a) 
(2.3b) 



where 4>l and (f>R are the chiral components of a Bose field : 4>l,b. — 1/2(0 =F / H(x')dx'). Thus there are two 
Bose fields <fi c and 0d corresponding to the two fermions ip c and ipd- The effective theory for the spin-1 chain can be 
simplified by use of the natural "acoustic" and "optic" linear combinations : 



The Bose theory can then be written as : 



where the a-sector is a sine- Gordon theory 
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and the o-sector is a generalized sine-Gordon theory 



Ho = 



K (Ul) + ^(V0 O ) 2 
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(2.4) 
(2.5) 
(2.6) 

(2.7) 



Here (f) — — </>jj Q is the dual field of <p , and the couplings are given by gi = g 2 = D — A, g 3 = —1 at first order in 
D, A. The velocities v a , v and the two parameters K a , K Q are also functions of the initial parameters of the problem 
and can be computed to the first order in D and A : 

3A + L> A — D „ 3A + D „ A-D 



1 
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Contrary to the case of the S=l/2 chain, there are no exact results on the values of the couplings K ai K a as functions 
of the bare parameters D and A. Strictly speaking, there are operators coupling the two fields 4> a and 4> in Eq.( [2.5| ) 
but they are less relevant than the interactions in Eq.(2.6) and Eq.(2.7). 

The physics of the a-sector is well-known : there is a Kosterlitz-Thouless phase transition separating a massless phase 
with algebraic quasi-long range order from a massive phase in which the field <j) a has a definite vacuum expectation 
value. The transition takes place when the operator cos(v / 87r0a) becomes relevant according to the value of the 
parameter K a : this happens for K a = 1. In the massive regime K a < 1 the field 4> a acquires a vacuum expectation 
value (4> a ) fixed by the minimum of the cosine operator appearing in the energy. In the massless regime, correlations 
functions decay algebraically with powers dictated by the value of the coupling K a . 

The physics of the o-sector is less transparent : it involves now the dual field <j) which is a non-local operator as 
a frmction of <p . In the Coulomb gas language, this theory has both electric and magnetic charges. In fact, den 
Nijsc.3 has shown that this effective theory appears also in the Coulomb gas representation of the XY model with Ising 
anisotropy. Since this model undergoes an Ising phase transition between two massive phases, we expect the same to 
be true for the generalized sine-Gordon theory ( ^.7[ ). 

Th e di fferent phases of the S=l spin chain can be characterized from the behavior of the effective theories (2.6) 
and (2.7). In the AF phase, both fields 4> a and 4> have a vacuum expectation value. In the Haldane and large-D 
phase it is now <j) a and <j) that have expectation values. It is the Ising transition of the theory (2.7) that takes place 



between AF and Haldane and also between AF and large-D. The difference between the Haldane and large-D phases 
is only the sign in front of the operator cos(v / 87T0a)- Thus it is the value of the condensate (<j) a ) that differs between 
the Haldane and large-D phases. When the coefficient of the operator cos(v / 87T0 a ) is zero, there is a massless line in 
the phase diagram separating H and D. When K a is large enough, the field 4> a is massless : when <fi condenses it the 
XY2 phase while 4> condenses in XY1 phase. These characteristics are summarized in table 1 (where "dis" means 
that the corresponding field has no mean value and is disordered). Finally, there is a ferromagnetic phase (F) where 
bosonization breaks down : at the F boundary the velocity v a vanishes. The F phase is beyond the scope of our study 
and will not be discussed. r-. 

The phase boundaries are in good agreement with the more de taile d picture of den Nijs and RommelseEI. It should 
be noted that the first-order expressions for K a and K a in Eqs. (|2.S| ) lead to straight lines for the phase boundaries 
while in reality there are tricritical points, i.e. between Haldane, large-D and AF phases. They are also out of reach 
of the bosonization method at the order which is considered here. 
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III. DISORDER IN THE BOSONIC LANGUAGE 



In this section we describe an approximate treatment of quenched randomness in the Bosonized theory obtained in 
sect. II. We give the operators that should be added to the Bosonic field theory ( |2.5| ) to describe randomness. We 
first list the various kinds of randomness that we have been able to treat. We will consider the effects of a random 
field along the z-axis. We add to the Hamiltonian a term : 

H ZF = Y,hiSl (3.1) 

i 

where the local fields hi are Gaussian random variables with zero mean, uncorrelated from site to site : 

hihj = VzF^j. (3.2) 

The bar means average over the random distribution and the variance of the distribution is T>zf, characterizing the 
strength of the disorder. This kind of disorder preserves the rotational symmetry around the z-axis which a symmetry 
of the spin hamiltonian. It is known from the Imry-Ma arguments that it has a dramatic effect on the phases with 
long-range magnetic order : an arbitrarily small amount of random z-field destroy the ordering. We consider also the 
case of a random planar field, i.e. lying in the XY plane : 

H PF = J2 h i S i +W, (3.3) 



with hfhj — T)pFSijS a p. The rotational symmetry is destroyed in this case. We treat random planar exchange : 



'Hpe — ^ Ji{SfSf +1 + S?S? +1 ), (3.4) 



with JiJj = VpESij. Finally randomness may break XY-symmetry by random planar anisotropy 



(3.5) 



with jijj = V PA 5 i: j. 

When bosonizing the spin operators in Eqs. ( |3~1 3~3],3.4,pT5|), there is one part with small momenta q w and one 
part with q « 7r. For example, in the ZF case we have : 



S z (x,t) 



cos(v 27t</) ) cos(v 2-K(j) ) 



(3.6) 



in the continuum limit. Inserted in Eq.(3.1) this formula will pick the q ~ part of the fluctuating field hi as well as 
the q ~ 7r part : 



Hzf= / dxh qa0 (x)O 1 {x,t) + h qa7T (x)O 2 (x,t), 



with 



0i = —\ —d x c/)a(x,t) and C 2 = cos(v 27r0 a ) cos(v 2n<j) ). 



(3.7) 



(3.8) 



In fact this formula is generic for all kinds of disorder. We have listed in Table II the operators 0\ and O2 corresponding 
to all cases cited above. We have kept only the most relevant operators in the renormalization group sense. The 
continuous random fields h q ~°(x) and /i 9 ~ 7r (x) are also Gaussian. To treat the quenched disorder we follow the replica 
trick : we introduce n copies of the system and compute the free energy as F = lim„^o {Z n — l)/«- Integrating over 
the Gaussian disorder h q ~°(x) leads in the Euclidean effective action to the following operator : 



S cff = -Vij dxdrdr' (x,T)o[ j) {x,t'), 



(3.9) 
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where T>i is the variance of h q ~°(x), i,j are replica indices. It is important to note that this operator is not local : 
since the random field h 9 (x) is quenched i.e. does not depend on time, the Gaussian integration leads to the 
two-time integral in Eq.( |3.9[ ). We obtain exactly similar formulas for all operators given in table II. In each case 
there is a variance T> corresponding to the Gaussian random field coupled to the operator considered. We will use 
the obvious notation ZF^°\ ZF^\ PF^°\ and so on. The disorder strengths are noted also T>z F , Dgp, .... The 
operator corresponds to a forward scattering process in the language of the Jordan- Wigner fermions. When 
written as a Bosonic operator it may be absorbed by a change of variable in the path integral over the field <f> a which 
leads to simple leading-order RG calculationsEj, when interactions in the a-sector are irrelevant. Here we will retain 
this Zf' ' operator explicitly in our RG treatment because, when its strength diverges, it leads to a very important 
phenomenon : the vanishing of the critical disorder strength for massive phases. 

It is important to note that the most relevant q = tt operators are, in all cases over study, always given by ZF^\ 
This is due to the fact that even when ZF^ does not appear at the bare level, it js generated by renormalization : 
this has been pointed out first by Doty and Fisher in the study of the S=l/2 chairJj. In our case a similar reasoning 
leads immediately to the appearance of ZF^ in the PF case as well as in the PA case (in the ZF and PE case, ZF^> 
appear at the bare level) 

The operators PA^ and PE^ are specific to the S—l spin chain as opposed to the spin ladder : they come from 
the cross coupling Oi ■ Ti+i that occurs when expressing S=l spin operators into S=l/2 entities. In the spin ladder it 
would be replaced by Oi ■ ti acro ss the rung. 

The operators of the type ( |3.9|) break the Lorentz invariance which is present in the effective field theory of the pure 
system. This means that in a renormalization group calculation many new operators will be generated in addition to 
those introduced at the bare level. In principle it is thus extremely difficult to keep the RG flow under control. In 
the RG calculation at second order that we performed, we have to add for internal RG consistency only three new 
operators : 

Sal = -Vai J dxdrdr' d T (P a {x,T)d T <jj a {x,T'), (3.10a) 

id 

Sa2 = -V a2 f dxdrdr' d w ^ a (x,r)^ a (x,r'), (3.10b) 

So = -V Q f dxdrdr' d z <p {x,r)ck4P {x,r r ), (3.10c) 
J i,3 

where z = x + v a r, z = x — v r, w = x + v a r and w — x — v a r. They are called al, a2 and o in what follows. 



IV. RENORMALIZATION GROUP TREATMENT 



In this section we perform a second-or der RG treatment of the Bosonic theory given by the field theory Eq. (2.5) 
perturbed by the operators of the type (|3.9|). There are five distinct bare operators (ZF^°\ ZF^\ PF^°\ PFW>, 
PA^) and three generated by renormalization (al, a2, o). The perturbation theory can start only from the massless 



theory with no cosine operators g\ = g% = <?3 = in Eq.(2.5). We write the effective theory including randomness as : 



S = S* +J2 a x '- 2 g t dxdrOi, (4.1) 



where S* is the massless theory characterized by K a and K a , a is the cut-off and Xi is the scaling dimension of the 
operator Oi . The beta functions of the couplings <?.; are then given by the following formula : 

^ = (2 - x k )g k - tt Y, C ijk g t 9j + 0(g 3 ), (4.2) 
where the Wilson coefficients Cijk arc obtained from the Operator Product Expansion : 

This procedure gives the beta function for an arbitrary number of replicas and is very useful in order to find which 
operators are generated by renormalization. While this technique can be applied straightforwardly to local operators 
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like those appearing in the sine-Gordon or generalized sine-Gordon theories, this is not so when dealing with the 
non-local operators ( |3.9|) generated by randomness. We will take into account only the strongest singularities (which 
are in l/\Rez\) in the OPE containing nonlocal operators, and assume that the OPE is still valid. 

We now sum up the renormalization group equations of the system. In the replica n — > limit, we obtain the flow 
for the cosine couplings : 

^ = 2(1 - K a ) Vl - -^K 2 ayi {2vf F + 2V al + V a2 ), (4.4a) 
^ = 2(1- K )y 2 -^K 2 V y 2 , (4.4b) 



where y\ = gi/v a , 2/2,3 = 92,3/ 'v . In the first equation (4.4a), in the absence of disorder the renormalization eigenvalue 



2(1 — K a ) leads to a massive phase when K a < 1 : the cosine is relevant since y\ scales to infinity in this case. This 
picture may be altered since the randomness leads to an additional second-order coupling —y\ x T^zf between the 
strength of the interaction and the forward scattering term ZF^°h If T>^ F (l) scales fast enough to infinity then it can 
revert the flow of y\ and leads to the destruction ofrthe massive phase. This behaviour has been obtained in the case 
of the Hubbard model by Fujimoto and KawakamiO. We will show later that this happens also in the case of the 
spin chain under study. Of course to capture such phenomena, it is important to deal explicitly with ZF^ instead 
of eliminating it by a change of variables in the functional integral. 
There are also two flow equations for the stiffness constants : 

dKa = -(2>g + ^vDK + \vf F + 2>& (4.5a) 



dl 



dK V M 2 X> (0) X> (0) 

ai\ n ' u ZF 1 * 2\ 7^2 1 i u PF 1 U PE , ± 2 



(^ + - 2 ylWi + {^ + ^ + - 2 yi), (4.5b) 

dl U 7T AU U 7T Z 

where u = v jv a . The operators that do not involve gradients have simple renormalization properties. In fact they 
scale according to the dimension that can be found immediately from the massless theory : 

d ^ = {i -K a -K )V^ F , (4.6a) 

^-(3-7^-7^)«, («b) 



dl v 4K a AK ' 



dV PE _ r , 1 y n - 



(4.6c) 



dl K K 



dV 



^ 1 

_|A = (3--pS?). (4.6d) 

These equations allow immediate statements about the stability of the massless phases according to the relevance or 
irrelevance of the random operators. This simplicity however does not extend to the operators ZF^°\ al, a2 and o 
because they appear in the OPE of various combinations of the other operators. The RG equations thus couple the 
various random perturbations : 

^ - ^ + ^\f^ V ZF 2 - ~ \hvf A \ (4.7a) 



G 



dV, 



1 

8?? 



1 

2u' 



(4.7b) 



Equations similar to (4.7) hold also for the couplings *D a \ aim 'Dai- The second-order terms in the flow equations 

n 

dy 



above involve the following functions 

dy 



/(«) 



(l + j / 2)K a /2( 1+u 2 y 2)K /2 



dy 



-i 2 



(1 +y 2)l/2K a 



fc(ti) 



4-oo 



( 1+y 2)l/8K a ( 1+u 2 y 2)l/8if 



dy 



{l + u 2 y 2 y/ 2K ° 



(4.8a) 



(4.8b) 



It is important to note that these functions are not well-defined for all values of the stiffness constants K a , K Q . We 
will show in the following section that these divergences lead to the disappearance of the Ising-ordered phase for an 
arbitrarily weak random z-field, a s predicted by the Imry-Ma argument. 

Armed with the flow equations (4.4,4.5 4. £,4. 7), we can study the phase diagram of the S=l chain in the presence of 
the various kinds of randomness listed in Table II. Before embarking in the general case, we briefly discuss a simplified 
model that display most of the behaviour of the complete theory : 

S eff = — I dxdT{d^4>f + y--^ / dx dr cos(V4ira 2 <j>(x,T))- (4.9) 



dx dr dr' cos(V 4irb 2 cj>(x, r)) cos(V Awb 2 (t)(x, t')). 



(4.10) 



Defining y = g/iru, and D — da/2itu 2 , where u is the velocity, we derive the following renormalization group equations 
for y, D and K : 



dy 

dl " 
dD 

~dT 

dK 
~df 



(2- 
= (3 



a 2 K)y, 
b 2 /K)D, 



-y 2 K 2 



D. 



(4.11a) 
(4.11b) 
(4.11c) 



The corresponding RG flow is given in Figure 2 when (ab) 2 < 6. Here, the crucial parameter is the stiffness constant 
K, because it governs the relative relevance of interaction and disorder operators. For instance, if K flows to zero for 
small initial disorder, the disorder operator will become irrelevant, whereas the interaction operator is relevant. If the 
pure case corresponds to a regime where interactions are relevant, the presence of infinitcsimally small (but finite) 
disorder will not affect this behavior. The value of K* in the pure case is enough to predict the qualitative behavior 
of the system under weak disorder at this order of perturbation theory. It is clear that the massive phase will persist 
up to some critical disorder which, strictly speakin g, is outside of reach of perturbation theory. A formal proof of 
stability is given in appendix A. The flow equations ( 4.11 ) are only first-order. The interest of our second-order study 
will appear in the study of the Ising-ordered phase. 

Finally, we comment on the correctness of the bosonization approach. Since the scheme we use is perturbation 
improved by use of the renormalization group then we can't say much about the flow to strong coupling. The phases 
that are unstable cannot be characterized by our tools. We can only find whether or not the randomness will have 
an immediate effect on the system. Even in the stable case, we of course expect that for strong enough disorder the 
physics of the system will change through a phase transition of some kind. This strategy has been applied with some 
success in the S=l/2 chain in a random z-field ans ip-ithe related problem of ID fermions with attractive interactions 
submitted to a random site potential. BosonizationE£tEI predicts stability of the XY phase,o£,the S=l/2 chain in the 
region —1 < A < —1/2 and there is ample evidence for the correctness of this resuliE3~E3 from direct numerical 
studies. Our general study should be followed by attacks by similar methods for an independent consistency check. 



V. PHASE DIAGRAM IN THE PRESENCE OF DISORDER 



In this section, we exploit the RG flow equations obtained previously to discuss the influence of randomness on the 
phase diagram of the S=l spin chain phase diagram pictured in Fig. (1). We begin our discussion by the massive 
phases, with and without long-range spin order. All our findings are summarized in Table III. The key element of the 
stability is the operator ZF^ since it is present in all kinds of randomness. 
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A. AF phase 



This phase is characterized by the following IR stiffness co nsta nts in the pure case : K* — K* = 0. For such a 
fixed point, the functions appearing at second-order in Eqs. (4.7) take the values g(u) = h = k(u) = while the 
function / (4.8a) is singular. The ZF^ operator is relevant in this phase. Since the cosine interaction is relevant in 
the a-sector, it seems that we are in a case with interaction and disorder that are relevant. However, the second-order 
beta functions allow a more precise conclusion : when K* and K* flow to a small enough value then / diverges. When 

approaching this singular point, T>^ F becomes arbitrarily large even for infinitesimally small values of V^ F {1 = 0) 
and T>^ F {1 = 0). In Eq. (4.4a), the coupling —y\ x V^ F becomes so large that the interaction y\ is driven to zero. 
We infer thus that the Mott-like gap of the Ising phase is immediately destroyed by disorder, i.e. by infinitesimal 
disorder. 

In the case of a random z-field, this is really what we predict from the Imry-Ma reasoning. In the case of the 
random planar held and random planar anisotropy, the random z-field is generated by the RG process so this is why 
there is also destruction of the Ising phase. The only trouble with the bosonization treatment is that it also leads to 
instability under random planar exchange, a fact which is not expected from the Imry-Ma argument. Since stability in 
this case is expected we see this as a problem of our approach (if we apply the same treatment to the S=l/2 case then 
one finds also instability under random planar exchange of the corresponding Ising phase). However this approach 
correctly reproduces the immediate destruction of the gap in the Ising phase. 



B. Haldane and large-D phases 



In the renormalization group equations, nothing distinguishes the H and D phases. So we expect that they behave in 
the same way under disorder. These phases are characterized by K* = 0, K* = oo, and <p a has a vacuum expectation 
value. The function / is now well defined when we reach the unperturbed fixed point, and near this fixed point, 
f(u) = 0. The operators ZF^\ PF^ ' and PA^ are irrelevant. This implies immediately stability under random 
XY symmetry-breaking exchange. Concerning the random fields, w e note that th e random forward scattering "D^f 
scales as e' while y\ and y 3 scale faster as e 2Z according to Eqs. (4. 7a, 4. 4a, 4.4c). Thus there is no possibility of 
destruction of the role of interactions contrary to the case of the AF phase. Thus there is stability under random 
z-field as well as under random XY-field. 

The only special case is the random XY-symmetric coupling. The PE^ operator is relevant and the function k(u) 
is singular. This drives the o-operator to large negative values of V a . This goes outside the reach of perturbation 
theory. But for such large negative values, y 3 scales much faster than the disorder like e^ D °^ 1 . We take this as an 
indication of the robustness of the gap against disorder. So we expect stability also in this case. The two gapped 
phases without long-range spin order are thus stable, up to some critical strength of the disorder presumably, which 
is beyond the reach of the methods we employ here. 



C. XY1 phase 

In this phase, the a-sector is gapless, and the ZF^ operator is irrelevant near the unperturbed fixed point in the 
whole phase so instabilities can appear only through the other operators induced by randomness. 

• Random z-field : the forward scattering ZF^ operator can be absorbed by a field redefinition. Thus the excitation 
spectrum is unchanged and the correlation functions are affected in a simple wayta. 

• Random XY field : the PF^ ' operator is relevant; so the phase is unstable. In this case, we note that the function 
g is singular, so the ZF^ operator is singular too and scale to +00. 

• Random XY symmetric coupling : the PE^ operator is relevant, as in the Haldane/large-D phase. The same 
discussion apply : the stability of the o-sector lead us to conclude that the XY1 phase remains stable under this kind 
of disorder. 

• Random XY symmetry-breaking coupling : The PA^ operator is relevant and the a-sector is not gapped; so the 
phase is unstable. The operator leading to the instability is different from the ladder case because it involves err 
couplings that are typical of the S=l case. 
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D. Frontier line between Haldane and large-D phases 

This line is characterized by yi — 0, so that the field <\> a is a free massless field. There is no gap along this line in 
the a-sector. It can be seen as an extension of the XY1 phase. Nevertheless, K a is renormalized by disorder and K a 
flows to infinity. Thus, the function g is equal to zero, whereas h is either finite or singular according to the value 
of K a . The stability discussion is exactly the same as that of the XY1 phase apart from the random XY field. The 
novelty is that K a may be less than 1 but the a-sector is still massless so there is a transition line in the T>p F — K a 
plane, corresponding to a critical value of (K a ) c =l/12. For K a > (K a ) c , the line we consider is unstable under weak 
PF disorder. The line starts from the XY1 boundary for which K a — 1 and crosses the frontier line between AF and 
H phases. The stiffness K a decreases along the line but we do not know if it reaches the value 1/12 before arriving at 
the tricritical point. We conclude that the line is unstable under a random XY field at least in the neighborhood of 
the XY phase. 



E. XY2 phase 



This phase is characterized by the following IR values in the pure case : K* > 1 and K* = 0. We are in a 
most interesting case in which the operator ZF^' may be either relevant or irrelevant within the bulk of this phase 



according to the scaling equation (4.6a). We first consider the influence of the random z-ficld : we just need to 
consider the a-sector and the operator ZF^\ which simplifies the calculation and does not change the result (the 
random forward scattering ZF^ does not change the global picture). The effective renormalization group equations 
are then at lowest order (since the function g is not singular near the IR pure fixed point, nothing new happens at 
the next order) : 

^ = 2(1 - K a )y u (5.1a) 



dV 



(7T) 

ZF _ (o _ u 



dl 



(3-K a )V%>, (5.1b) 



d -^r--^l + i vl)Kl (5 - lc) 

In each plane y\ — and T) ZF — 0, the flow has a simple form. We can draw the corresponding renormalization flow 
in the three dimensional space : it is given in Figure 3. We observe that the ZF^ operator breaks the XY2 phase 
into two phases. The first one (for large K) is a massless phase stable under small disorder, and the interaction term 
also flows to zero; the other phase is unstable under weak disorder. Thus the XY2 phase is only partially stable under 
small disorder, and this is consistent with physical intuition. Indeed, when the single-ion D(S'f) 2 term becomes large 
and negative, we expect that the S=l spin chain will behave as the S=l/2 spin chain. As shown by Doty and Fisher 
for the S=l/2 spin chain, there is a region of exchange parameter A stable under weak disorder (it is also a superfluid 
phase arising for the disordered boson gas; see Giamarchi and Schulz). In the S=l/2 spin chain, this stability region 
is located in the interval — 1 < A < —1/2. In the S=l case, we find that such a stable phase also exists for K large 
enough, and this is consistent with the fact that we should recover the S=l/2 behavior. 

In the case of a random XY field, since there is a random z-field generated we are in the same situation as above. 
For the random XY symmetric coupling, the PE^ operator is irrelevant and the above discussion is again valid. 

In the random XY symmetry-breaking case, the PA^ operator is always relevant whereas the a-sector is ungapped. 
Thus, the chain is unstable. This is again consistent in the limit of large negative D with the Doty and Fisher results 
for the S=l/2 chain. 

Finally, we briefly comment the case of random z-exchange that we have not considered yet in this paper. In terms 
of the S=l spins, it is given by : TtzE = J2i Ji^i^i+v When written with the two kinds of spins S=l/2, terms of the 
form err appear. They lead to cos(\/87r0a) + cos(\/87r0o) when expressed in the boson language, a coupling that has 
not been studied in the context of spin ladders. The cos(-\/87r0a) leads to an operator which is relevant in Haldane 
phase, and make the forward scattering ZF^ singular and infinite with the consequence of vanishing interactions in 
the a-sector. The AF phase is also unstable : the ZF^ operator is also contained in the boson expression of this 
kind of disorder, as if a randora, z-field was generated. These instabilities clearly deserve more studies. This problem 
appears also in the S=l/2 cases. The XY1 phase is stable while XY2 is unstable. 
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VI. CONCLUSION 



We have performed a renormalization group study of the influence of disorder on the effective theory which describes 
the S=l chain. The RG equations allow a discussion of the stability of the various phases of the system. Our discussion 
includes the effects of the second order of the renormalization group calculation. This allows us to capture the Imry-Ma 
destruction of the AF phase. 

The gapped Haldane and large-D phases are stable under all kinds of disorder we have studied: random z-field, 
random XY field and random XY symmetric and antisymmetric exchange. On the contrary, the AF phase which is 
also gapped is less stable. In fact, our RG calculation that goes up to the second order in the disorder strength shows 
that the phase is unstable to an arbitrarily weak random z-field, in agreement with the Imry-Ma argument. It is 
also unstable under random XY-field and random planar anisotropy since these perturbations do generate a random 
z-field. 

Concerning the gapless XY phases, then again the situation is quite rich. The XY1 phase is stable under random z- 
field and random planar exchange while it is immediately unstable under perturbations breaking the planar symmetry : 
random XY-field and planar anisotropy. The XY2 phase is only partially stable under random fields and random 
planar exchange. In this case the phase breaks up into two parts: there is a stable massless domain with irrelevant 
disorder and another domain with relevant disorder. Finally XY2 is totally unstable under random planar anisotropy. 

With respect to the closely related problem of 2-leg spin ladder, the effective theory is the same but the the PE 
and PA type of disorder leads to different operators in the bosonization approach. As a consequence, there are some 
instabilities that are induced by a different mechanism : this happens in the massless phases XY1 and XY2 and also 
along the line separating the Haldane and large-D phases. 

Of course it remains to characterize more completely the phases in which the disorder is relevant. In our calculation, 
we can simply observe that the system flows to strong coupling but the methods we use are not informative on its fate. 
An outstanding problem is to describe in a unified manner the weak-coupling regime we observe here with the regime 
probedLJay the real-space RG calculations. The gapless phase with hidden long-range order of the random spin-1 
chainEjIlZI appears presumably only beyond some critical strength and it seems to be out of reach of the methods we 
have used in this paper. 

Note added: A recent work by Y. Nishiyama on the effect of a random z-field on the Haldane gap and the XY1 



phase is in agreement with our findings (see e-print cond-mat/9805110) 
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APPENDIX A: STABILITY OF A GAPPED PHASE 

The pure SG system is gapped when K — ► 0. We call K° and y° a solution of the pure differential system. We 
set K = K° + SK and y = y° + Sy. A straightforward calculation leads to : 

SK' = D — 2y a (K°) 2 Sy - 2(y°) 2 K°SK (Al) 
Sy' = (2 - a 2 K n )Sy - a 2 y°SK (A2) 

From the flow equations, we have D(l) < Doe 31 , SK > and Sy > 0. Furthermore, the initial values of the functions 
SK and 5y are 0. From this statement, we can solve exactly the differential system for SK and Sy if we replace D(l) 
by D e 31 . We get : 

SK(l) = D F(l) (A3) 
Sy(l) = D G(l) (A4) 

where the functions F and G obey the following differential system : 

F'(x) = e 3x - 2y (x) 2 K°(x)F{x) - 2y°(x)K°(x) 2 (A5) 
G\x) = (2 - a 2 K°(x))G(x) - a 2 y°(x)F(x) (A6) 
F(0) = G(0) = 0. (A7) 
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If we choose a small e > 0, it is possible to find I = Iq so that K°(lo) < e/2. Moreover, we know that the real 
function 8K is inferior to DqF(x) where F(x) is a known function; so we can choose the initial value of D(l) so that 
5K(Iq) < e/2. In this case we have K(Iq) < e. Because of the positivity of the functions SK and Sy, we know that 
y 2 (l)K 2 (l) > (y°) 2 (l)(K Q ) 2 (l) for all 1, and so dK/dl < ~(y°) 2 (l)(K°) 2 (l) + D e 31 . This inequality proves that for 
sufficiently small Do, we can also have K'(Iq) < 0. This argument shows that, for all e > 0, there is a positive Iq and 
a small enough D so that K (Iq) < e and K'Qq) < 0. For / = Iq, 3 — b 2 /K(l n ) ~ —b 2 /K(lo). As a result, the disorder 
contribution in K' will decrease to as fast as e~ b / K W l , whereas the other term will decrease slower than K 2 (l). 
Thus this last term dominates the value of K' which remains negative for all I > Iq. It follows that for I 3> Iq, we 
have K' ~ —y 2 K 2 , so that K(l) — > for a sufficiently small but finite value of Dq. 
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FIG. 3. Effective renormalization flow for AF and XY2 phases 
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TABLE I. Values of the different parameters characterizing the phases of the pure S=l antiferromagnetic chain. 





q = 


q — ir 


ZF 


d x (f>a 


cos(\/2ir(f) a ) cos(\/2n(f>o) 


PF 


COs(y / n/2<j> a ) COs(y / V /2(j> ) 


cos(\/2~ir(f)a) cos(^/27r0 o ) 


PE 


COS(V27T0 O ) 


cos(\/27r0 a ) cos(\/27r0 o ) 


PA 


cos(\/27r0a) 


cos(\/27r0 a ) cos(V27T(?!) ) 



TABLE II. Bosonized formulas of the quenched disorder. 





AF 


Haldane 


large- D 


XY1 


XY2 


ZF 


Unstable 


Stable 


Stable 


Stable 


Partially stable 


PF 


Unstable 


Stable 


Stable 


Unstable 


Partially stable 


PE 


Unstable 


Stable 


Stable 


Stable 


Partially stable 


PA 


Unstable 


Stable 


Stable 


Unstable 


Unstable 



TABLE III. Stability of the S=l phases with respect to quenched disorder. 
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